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I ' Abstract 

We consider the possibility of measuring both sin(2/3) and cos(2/3) in the KM unitarity triangle using 
the process B — > D* + D*~ K s . This decay mode has a higher branching fraction (0(1%)) than the mode 
£>° — > D* + D*~ . We use the factorization assumption and heavy hadron chiral perturbation theory to 
estimate the branching fraction and polarization. The time dependent rate for B°(t) — > D* + D*~ K s can 
be used to measure sin(2/?) and cos(2/3) . Furthermore, examination of the D* + K s mass spectrum may 
be the best way to experimentally find the broad 1 + p-wave D s meson. 



1 Introduction 



The decay B° — ► J /ipK s is expected to provide a clean measurement of the angle sin(2/3) in the unitarity 
triangle Q. However, other modes can also provide relevant information on the angle (3. An example of 
such a mode is the decay B° — > D^D . In this mode -B — ► D* + D*~ , the vector- vector final state in 
general is an admixture of CP odd and even eigenstates, because s, p and d partial waves with different 
CP-parities can contribute. Since the CP asymmetry has opposite sign for the two CP states they tend 
to cancel or dilute the overall asymmetry. The amount of dilution of the CP asymmetry is represented 
by the dilution factor, D, which depends on the CP composition of the final state. The presence of two 
CP components in the final state of B° — > D* + D*~ makes the dilution factor, D < 1 for this decay. This 
is unlike the case for a mode such as B° — > D + D~ where the final state is a CP eigenstate and D=l as 
there is no dilution of the CP asymmetry. An angular analysis can extract the contribution of the different 
CP eigenstates leading to a measurement of D and hence of sin(2/3) [Q, B. However, in the factorization 
approximation and using Heavy Quark Effective Theory (HQET) it can be shown that the final state in 
B° — > D* + D*~ is dominated by a single CP eigenstate @. To the extent that this is valid, the angle 
sin(2/3) can be determined without the need for an angular analysis. The decay B° — > D* + D*~ may be 
preferred to B° — > D + D~ because contamination from penguin contributions and final state interactions 
(FSI) are expected to be smaller in the former decay ||. 

In this work we consider the possibility of extracting j3 from the decay B° — ► D(*'L)(*'K S . These modes 
are enhanced relative to B° — >■ D^D^ by the factor \V cs /V C( i\ 2 ~ 20. As in the case of B° — > J/i>K s 
decay, the penguin contamination is expected to be small in these decays. Moreover these decays can be 
used to probe both sin 2/3 and cos 2/3 which can resolve (3 — ► ir/2 — (3 ambiguity ||. 

The possibility that a large portion of the b — ► ccs decays materialize as B — > DDK was first suggested 
by Buchalla et al. Q. Using wrong sign D-lepton correlations, experimental evidence for this possibility 
was found by CLEO, who observed B(B -> DX) = (7.9 ± 2.2)%Q. Later, CLEO§, ALEPH|] and 
DELPHl|lO| reported full reconstruction of exclusive DDK final states with branching fractions that are 
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consistent with the result from D-lepton correlations. 

CLEO obtained B(B° -> D* + D* R-) = (1.30t£g ± 0.27)% and B{B~ -► D* D* K~) = (1.45^ ± 
0.36)%. These values should be approximately equal to the branching fraction for B(B° — > D* + D*~ K°). 
We use the latter value for the purpose of a sensitivity estimate. Taking into account B(K° — ► ET S ) = 0.5, 
B(K S — > 7T + 7r _ ) = 0.667, and assuming that the if s reconstruction efficiency is ~ 0.5, we can estimate the 
ratio of the tagged B° — > D* + D*~ K s events to the tagged D* + D*~ events. Assuming — > D* + D*~) = 



6 x 10~ 4 , which is the central value of the recent CLEO measurement [11], we find that the ratio of the 
number of events is ~ 4.0. Therefore, this mode could be more sensitive to the CP violation angle sin(2/3) 
than B° —* D* + D*~ . However, if the final state contains a resonance then B° and B° can be distinguished 
and there is additional dilution of the CP asymmetry. For the decay B — > f and B — > f the dilution 
factor, D, measures the ratio of the overlap of the amplitudes for B — > f and B — ► / to the average of the 
decay rate for B —> f and B — ► / . Clearly D = 1 when the amplitudes for B — > f and B — ► / decays are 
equal. When the final state in the decay B — ► D* + D*~ K s contains a resonance the amplitude for B and -B 
decays are different because the resonance in the B and B final state occur at different kinematical points. 
This causes additional mismatch of the B and B amplitudes which results in the further dilution of the 
CP asymmetry. A similar conclusion is obtained in the comparison of B° — > D + D~K S to B° — > D + D~ . 
The above conclusions are detector dependent; a somewhat pessimistic estimate of the -fT s reconstruction 
efficiency is used here while the detection efficiency for the D* + D*~ final state is assumed to be similar 
for both cases. Better determination of the CP sensitivities will require more precise measurements of 
the branching fractions for the D*D*K decay modes and will also depend on details of the experimental 
apparatus and reconstruction programs. 

The amplitude for the decay B° — > D*D*K S can have a resonant contribution and a non-resonant 
contribution. For the resonant contribution the D*K S in the final state comes dominantly from an excited 
D s (l + ) state. In the approximation of treating D*D*K S as D* D s (ex cited), there are four possible excited 
p-wave D s states which might contribute. These are the two states with the light degrees of freedom in a 
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j p = 3/2 + state and the two states with light degrees of freedom in a j p = l/2 + state. Since the states 
with j p = 3/2 + decay via d-wave to D*K S , they are suppressed. Of the states with light degrees of freedom 
in j p = l/2 + states, only the 1 + state contributes. The + state is forbidden to decay to the final state 
D*K S . 

To estimate the above contribution and to calculate the non-resonant amplitude, we use heavy hadron 
chiral perturbation theory (HHCHPT) [|i~2"| . The momentum of K s can have a maximum value of about 



1 GeV for B° — > D* + D*~ K s . This is of the same order as A x which sets the scale below which we expect 
HHCHPT to be valid. It follows that in the present case it is reasonable to apply HHCHPT to calculate 
the three body decays. 

In the lowest order in the HHCHPT expansion, contributions to the decay amplitude come from the 
contact interaction terms and the pole diagrams which give rise to the non-resonant and resonant contri- 
butions respectively. The pole diagrams get contributions from the various multiplets involving D s type 
resonances as mentioned above. In the framework of HHCHPT, the ground state heavy meson has the 
light degrees of freedom in a spin-parity state j p = \ , corresponding to the usual pseudoscalar-vector 
meson doublet with J p = (CP, 1~). The first excited state involves a p-wave excitation, in which the light 
degrees of freedom have j p = ^ + or | + . In the latter case we have a heavy doublet with J p = (l + ,2 + ). 
These states can probably be identified with .D s i(2536) and D s j(2573) fl~3]| . Heavy quark symmetry rules 
out any pseudoscalar coupling of this doublet to the ground state at lowest order in the chiral expansion 



[14 1; hence the effects of these states will be suppressed and we will ignore them in our analysis. In fact 
there is an experimental upper limit on inclusive B -> L> s i(2536)X < 0.95% at 90% C.L gjj. Since the 
total D*D K rate is about 8%, this confirms that the narrow p-wave states do not account for a significant 
fraction of the total D*D*K rate. 

The other excited doublet has J p = (0 + ,l + ). These states are expected to decay rapidly through 
s-wave pion emission and have large widths |i~q| . Observation of the 1 + state in the D system was recently 



reported by CLEO 17]. Only the 1 + can contribute in this case. For later reference, we denote this state 
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by D*'i- However, quark model estimates suggest |l8|] that these states should have masses near m + dm 
with 5m = 500 MeV, where m is the mass of the lowest multiplet. 

We will assume that the leading order terms in HHCHPT give the dominant contribution to the decay 
amplitude and so we will neglect all sub-leading effects suppressed by 1/A X and 1/m, where m is the heavy 
quark mass. We show that from the time dependent analysis of B°(t) — ► D* + D*~K S one can extract 
sin(2/3) and cos(2/3). Measurement of both sin(2/3) and cos(2/3) can resolve the (3 — > ir/2 — (3 ambiguity 
]|, Il9|, [2(1. The measurement of sin(2/3) can be made from the time dependent partial rate asymmetry 
while a fit to the time dependent rate for T[B°(t) -> D*+.D*~.fir s ] + T[B°(t) -> £)*+£)*- i^ s ] may be used 
for the extraction of cos(2/3). Note that the cos (2/3) term measures the overlap of the imaginary part of the 
amplitudes for B — * D* + D*~ K s and B — * D* + D*~ K s decays and is non zero only if there is a resonance 
contribution. 

As in the case for B — > D* + D*~ the asymmetry in B —* D* + D*~ K s is also diluted. For the non 
resonant contribution to B — > D* + D*~ K s the final state is an admixture of CP states with different CP 
parities. This leads to the dilution of the asymmetry and this is the same dilution of the asymmetry as 
in the case for B — > D* + D*~. As already mentioned above, when the resonant contribution is included 
there is further dilution of the asymmetry from the additional mismatch of the amplitudes for B and B 
decays. One can reduce the additional dilution of the CP asymmetry by imposing cuts to remove the 
resonance. A narrow resonance is preferable as it can be more effectively removed from the signal region 
than a broad resonance . In this work we examine several cuts that can be used to remove the resonance 
and lessen the dilution of the CP asymmetry. When we include the resonance contribution we find that 
a broader resonance leads to a larger value of D and is a more useful probe of cos (2/3) because of the the 
larger overlap of the amplitudes for B — > D* + D*~ K s and B — > D* + D*~ K s decays. 

We also point out that from the differential decay distribution of the time independent process B° — ► 
D* + D*~ K s one can discover the 1 + resonance D*[. We show that the differential decay distribution for 
small values of E^, the kaon energy, shows a clear resonant structure which comes from the pole contribution 
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to the amplitude with the excited J p = 1 + intermediate state. Therefore, examination of the D*K S mass 
spectrum may be the best experimental way to find the broad 1 + p-wave D s meson and a fit to the decay 
distribution will measure its mass and the coupling. 

A similar analysis can be performed for B° — > D + D~ K s M, 21]. However the predictions of HHCHPT 



for this mode may be less reliable because of the larger energy of the K s . The effects of penguin contribu- 
tions, though small, may also be more important in B° — ► D + D~ K s than in B° — > D*D*K S as in the two 
body case || . 

In the next section we describe the extraction of of sin 2(3 and cos 2(3 from the time dependent rate for 
B{t) — ► D* + D*~ K s . In the next section we present the the amplitude for B — > D*D*K S in the factorization 
approximation and using HHCHPT. In the last section we discuss and present our results. 

2 Extraction of sin 2/3 and cos 2/3 

In this section we discuss the extraction of sin 2(3 and cos 2(3 from the time dependent rate for B(t) — > 
D* + D*~ K s . We define the following amplitudes 



a 



XiM = A[B o {p) _ D +*{p + )D^(p^K s (p k )), a x ^ = A(B\ V ) - D^(p + )D^(p^K a (p k ), (1) 



where B° and B° represent unmixed neutral B and Ai and A2 are the polarization indices of the D* + and 
D*~ respectively. 

The time-dependent amplitudes for an oscillating state B°(t) which has been tagged as a B° meson at 
time t = is given by, 

A x ^(t) = a ^ cos (^l)+ie- 2 ^a x ^ S m(^-) , (2) 



and the time-dependent amplitude squared summed over polarizations and integrated over the phase space 
angles is: 

\A(s + ,s~;t)\ 2 = - [G (s + , ) + G c (s + , s" ) cos Am ( - G s ) sin Am t] (3) 



with 

G (s + ,O = |a( S + ,0| 2 + |a( S +,0| 2 , (4) 

G c ( S +,0 = \a(s + ,s-)\ 2 -\a(s + ,s~)\ 2 , (5) 

G s (s + ,s-) = 29(e- 2i/3 a(s+s-)a*(s + ,s-)) 

= -2 sin(2/3) 5ft (aa*) + 2 cos(2/3) 9 (aa*) . (6) 

The variables s + and s~ are the Dalitz plot variable 

s + = (p++Pk) 2 , s~ = (p-+p k ) 2 

The transformation defining the CP-conjugate channel §°(t) — > D*~ D* + K s is s + <-> s _ , a <-> a and 
/3 -► -/3. Then: 

|A(s",s + ;t)| 2 = - [G (s",s + ) - G c (s",s + )cos Amt + G s (s~, s + ) sin Ami] . (7) 

Note that for simplicity the e~ rt and constant phase space factors have been omitted in the above equations. 

It is convenient in our case to replace the variables s + and s~ by the variables y and where E k is 

the K s energy in the rest frame of the B and y = cos 9 with 9 being the angle between the momentum of 

K s and D* + in a frame where the two D* are moving back to back along the z- axis. This frame is boosted 

with respect to the rest frame of the B with (3 = — (pfc/ms)(l/(l — Ek/ms))- Note s + <-> s~ corresponds 

to y <-> — y. The variable y can be expressed in terms of variables in the rest frame of B. For instance 

E' B E' + -p' B p' + y 

h + = 

vriB 

where E + and E' + are the energy of the D* + in the rest frame of the B and in the boosted frame while 
E' B is the energy of the B in the boosted frame. The magnitudes of the momentum of the B and the D* + 
in the boosted frame are given by p' B and p' + respectively. 

In the approximation of neglecting the penguin contributions to the amplitude there is no direct CP 
violation. This leads to the relation 

a x ^(p kl ,E k ) = a- x ^- x *(-p kl ,E k ) (8) 



where p k \ is the momentum of the of the K s in the boosted frame. The above relations then leads to 

G (-y,E k ) = G (y,E k ) (9) 

G c (-y,E k ) = -G c (y,E k ) (10) 

G s i(-y,E k ) = G sl (y,E k ) (11) 

G s2 (-y,E k ) = -G s2 (y,E k ) (12) 

where we have defined 

G sl (y,E k ) = »(aa*) (13) 

G s2 (-y,E k ) = Q(aa*) (14) 

Carrying out the integration over the phase space variables y and E k one gets the following expressions for 
the time-dependent total rates for B°(t) — > D* + D*~ K s and the CP conjugate process 

r(t) = -[I + 2sin(2/3)sin(Amt)/ s i] (15) 
r(t) = i[/ -2sin(2/3)sin(Amt)/ sl ] (16) 

where Jo and I s i are the integrated Go(y, E k ) and G s \(y, E k ) functions. One can then extract sin(2/3) from 
the rate asymmetry 

r(*) - r(t) _ D gin(2 ^ sin(Amt) (17) 



T(t) + T(t) 



where 



is the dilution factor. 

The cos (2/3) term can be probed by by integrating over half the range of the variable y which can be 
taken for instance to be y > 0. In this case we have 

r(t) = 7^0 + J c cos(Amt) + 2sin(2/3)sin(Amt)J s i - 2 cos(2/3) sin(Amt) J s2 ] (19) 
T(t) = ^[J + J c cos(Amt) - 2sin(2/3)sin(Amt)J s i - 2cos(2/3) sin(Ami) J s2 ] (20) 



where J , J c , J s i and J s2 , are the integrated G (y,E k ), G c (y,E k ), G s i(y,E k ) and G s2 (y,E k ) functions 
integrated over the range y > 0. One can measure cos(2/3) by fitting to the time distribution of T(t) + T(t). 
Measurement of the cos(2/3) can resolve the (3 — ► j — (5 ambiguity. 

3 Amplitude and Decay Distribution 

In this section we present the amplitude and decay distribution for the decay B —* D* + D*~ K s . Details of 
the calculation of the amplitudes using the factorization assumption and HHCHPT are given in Appendix 
A. 

The non-resonant amplitude for the three body decay B°(v, m) — > D* + {e\, v + , mi)D*~ (e 2 , u_, m2)K s {p k ), 
after setting m.2 = mi, is given by 

Id* 



Inon-res = K y/m^/m[mi£{v ■ V+) 1 — 

JK 

is^el^vav+p + e\ ■ ve* 2 • v+ - e\ ■ e* 2 (v • v + + 1)] (21) 

where K = + C2^j ■ Note that the amplitude above is the same as the amplitude for B° — > D*~ D* + 

Q except for a constant multiplicative factor ~ 1//at- 

To a good approximation one can use v ~ where v is the velocity of the B° in the boosted frame 
where the two D* are moving back to back. The K s , in this limit, is emitted in a s-wave configuration as 
the amplitude is independent of the angles that specify the K s momentum in the boosted frame. Then, 
as in the B° — > D* + D*~ case there are three helicity states allowed, (+,+), (— , — ) and (0,0) with the 

corresponding helicity amplitudes -£/++, H and Hqq. The helicity states are not CP eigenstates but one 

can go to the partial wave basis or the transverse basis where the states are CP eigenstates. The transverse 
basis amplitudes are related to the helicity amplitudes as 



A 
A 



V2 
H4-1- - H. 



V2 

A = H 00 (22) 



The three partial waves that are allowed in this case, s, p and d are then given by 

V2A» - A 



s = 



d 



V3 

p = A ± 
V2A + A\\ 

7T~ 



(23) 



The CP of the final state is given by ri(—) L where rj is the intrinsic parity of the final states and L is 
the relative angular momentum between D* + and D*~ . In the approximation v ~ one can write the 
non-resonant amplitude for B°(v,m) — > D* + (ei, v + , mi)D*~ (e 2 , V-, mi)K s (pj c ) 



'■n on— res 



JK 

-ie^e^e^vav-p + e* 2 ■ ve\ ■ v- - e[ ■ e* 2 (v ■ v- + 1) 



(24) 



There can also be pole contributions of the type shown in Figure 1. 



B 




Figure 1: The pole contribution to the process B — > D*D*K$- The intermediate state / can be D*^ or 
D*. The solid square represents the weak vertex while the solid circle represents the strong vertex. 

These give the decay sequences 



and 



B° -> D* + D*'{~ -> D* + D*~K° 



B° -> D* + D*~ -> D* + D*~K° 
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The propagator for the vector resonance is given by 

_ KVpVv - g^ u ) 



(25) 



2V ■ k 

where the momentum of the propagating particle P = m{V + k where to/ is the mass of the intermediate 
particle in Figure. 1. 

The contribution from the pole diagrams are given by ai res and ~a~2res, where a\ res is, with mi = m*' , 

fD*[ hp k ■ V- 



K\fm^/m\^/m\ v m*'£(v.v + ) 



(pfc • w_ + mi — m*' H 2 sl - 



ie^e\^* 2u v a v +p - e\ ■ ve* 2 ■ v + + ej ■ e* 2 (v ■ v + + 1). 



(26) 



Note that the above amplitude can be rewritten as 



Id*[ /to* 



hpk • v. 



*non—res 



V mi . v _ _|_ mi _ m *l _|_ 

a 2res is given by, with mi = m* where m* is the 1~ D* mass, 

. Id% g 



si 



(27) 



«2r 



Kyfm y Jmi^miyJm*£ i {v.v + ) 



f R (p k ■ V- + (mi - to,*) + 



-X 



X = -ie^ val3 e* 2ll p ku v +a v_p€\.v + ie^^el^p^v^v ■ v + + 1) + (e* • v_e2 • upfc ■ v+ - e* ■ v„e* 2 ■ v + p k ■ v) 
+ (e* • Pfc£2 • u + w •«--€*■ p fc €2 • w+ • «-) + (e* • e* 2 p k • ot+ • u_ - e£ ■ e* 2 p k ■ v + v • v_) (28) 

The amplitude <22res gives a tiny contribution to the total amplitude and can be neglected. In fact, 



this amplitude vanishes in the small velocity limit where the D* are almost at rest [22]. We note that the 
process with the + intermediate state 



B° -► D* + D; -> D* + D*-R° 

is not allowed due to parity conservation while the amplitude with the 0~ intermediate state 

B° _> D*+DJ -> D* + D*-R° 

is expected to be small compared to ai res . The propagator term in the above amplitude goes as approxi- 
mately \/{Ek + (to_d* — rriD a )) which does not have a pole as in aires- Moreover, the amplitude is further 
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suppressed with respect to ai res by a factor ~ Pk/Ek or |w|/t>o, where v and vo are the three velocity and 
the time component of the velocity four vector of the D* , from the DfD* + K° vertex. 

The total amplitude for B°(v,m) — > D* + (ei,v+,mi)D*~ (e2,v-,mi)K s (p k ) can be written as 

a = a non - res [l - Pi] (29) 

and the total amplitude for B°(v,m) — ► D* + (ti,v + ,mi)D*~ (e2,v-,mi)K s (p k ) can be written as 

a = a non - res [l - P 2 ] (30) 

with 



f D *[ m*' hp k ■ v 

1 ~ 777 V 777 (31) 

p 2 = ^ hp k ■ v+ ^ 

Id* V m i ( pfe . v+ + mi _ m */ + 

Note Pi and P2 can be expressed in terms of Pfc and y and P\{y,Ek) = P2(—y,E k ). The relation 
between quantities in the boosted frame and the rest frame of the B along with the calculation of the 
squared amplitude are given in Appendix B. 

The double differential decay distribution for the time independent process 

B°(v,m) -»■ D* + (ei,v + ,m 1 )D*-(e 2 ,V-,m 1 )K s (p k ) 

can be written as 

1 dT = f(y,E k ) 
TdydE k j f{y , Ek f_k± dydEk 

where p' k and p' + are the magnitudes of the three momentum of the K s and D* + in the boosted frame and 

f *' 

the expression for f(y, E k ) can be found in Appendix B. The differential distribution depends only on -j^, 
the mass m*' and the coupling h of the D*[ state . It is expected that /£>*/ « /d* and in the 5C/(3) limit 
Id* = Id* ■ So in the SU (3) limit a two parameter fit to the differential decay distribution can determine 
the mass and the coupling of the D*[ state . 
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The widths of the positive parity excited states are expected to be saturated by single kaon transitions 
14f. In our calculation we require the width of the D*[ state. Assuming 



T D * +/ « T(D*J' -> D* + K u ) + T(D*J -» D W K + ) (34) 

si 



one can write 



r D ., = -A^-I^m*' -mi) 2 p (35) 
tt/^- to*' 



where p is the magnitude of the three momentum of the decay products in the rest frame of D*[ and mi 
and m*' are the masses of the D* and D*[ state . 

It is clear that if a = a then the dilution factor D = 1. However that is not the case here. For 
the non resonant contribution, in the approximation of small velocity of the B, the final state is an 
admixture of CP states with different CP parities. This leads to D < 1. This is the same dilution of the 
asymmetry as in the case for B — > D* + D*~ . When the resonant contribution is included the amplitude 
a and a have an asymmetric dependence on the variable y. This reflects the fact that in the process 
B° — » D* + D*'{~ — > D* + D*~K° the kaon emerges most of the time closer to D*~ than the D* + . The 
situation is reversed for B° decays. Consequently there is additional mismatch between the amplitudes a 
and a which leads to further dilution of the asymmetry. One can reduce the dilution of the asymmetry, 
i.e., increase D, by imposing cuts so as to reduce the resonant contribution. We consider several cases 
where cuts may be employed to decrease the dilution of the asymmetry. From Eq. (31-32) it is clear that 
resonance occurs when the following condition is met 

Pk • v+ = m* — mi (36) 
Pk ■ v_ = m* — mi (37) 

If, in the allowed region of Ek, we can find a value Eko such that for values of E^ > Eko the above 
conditions are not satisfied for — 1 < y < 1 then we can remove the resonance by using the cut E^ > E^q . 
The value of E^q ~ 0.76 GeV in our case. We will call this case cut 1 for future reference. 
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Another possible cut is to include the whole range of but in the region E\. < E^q we remove the 
resonance by cutting on the variable y. We can use the region —0.5 < y < 0.5 since for most values of 
Ek the resonance condition is satisfied in the range — 1 < y < —0.5 and 0.5 < y < 1. We will call this 
case cut 2 for future reference. In any event, the cuts can be optimized after the resonance has been seen 
experimentally. However as we try to increase the value of D by reducing the resonance through cuts we 
also lessen the branching ratio from loss of signal. 

4 Results and Discussions 
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Figure 2: The branching fraction for B° — > D* + D*~ K s as a function of the h with and without cuts. 

As inputs to the calculation, we use fo% ~ = 200 MeV and take the mass of the D*[ state to be 
2.6 GeV. For the Isgur-Wise function we use the form 



QCD sum rules have been used to compute the strong coupling constants g and h |2jJ. We will use g = 0.3 
as obtained in Ref [^] but keep h as a free parameter because this coupling plays a more important role 
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Figure 3: The Dilution factor D as a function of the h with and without cuts, 
in the decay widths. 

Fig. 2 shows the branching fraction for B° — » D* + D*~ K s as a function of the coupling h. A QCD 
sum rule calculation gives h ~ —0.5 |p3[ |. We use the same sign of h as obtained in QCD sum rule 
calculation but vary h from —0.6 to —0.1. For this range of h the branching fraction can vary in the range 
0.45 — 0.93% when we employ no cuts. For h = —0.4 which corresponds to a D*[ state with a width 
of about 150 MeV the branching fraction is 0.83%. In our calculation this corresponds to a branching 
ratio B(B° -> D*-D*+K°) B(B° -> D*-D*°K+) « -» D*°D*°K + ) « £(5+ -> D*°D*+K°) 

0.9 — 1.86%. This is consistent with the CLEO measurements mentioned above. Fig. 2 also shows the 
branching fractions after cuts have been applied to reduce the dilution of the CP asymmetry. 

Fig. 3 shows a plot of the dilution factor D versus the coupling h. In the absence of any cuts we find 
that larger values of \h\ gives a larger value of D and hence less dilution in the asymmetry because for a 
broad D*[ state there is more overlap between the amplitudes for B° — ► D* + D*~ K s and 1?° — > D* + D*~ K s . 
For /i = —0.4 the dilution factor is about 0.75 with no cuts. For the case of cut 1, where we use the cut 
Ek > E^o to effectively remove the resonance, the dilution factor increases with smaller \h\. This is because 
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Figure 4: The squared amplitude for B° — > D* + D* K s B° — > D* + D* K s as a function of the variable y 
for /i = —0.4 which corresponds to a D*^ state with a width of about 150 MeV. 

for smaller \h\ and Ej~ > E^q the resonant amplitude is small and the total amplitude is dominated by the 
non-resonant amplitude which give a larger value for D. For the case of cut 2, as in the case with no cuts, 
the dilution factor D decreases with smaller \h\. This is because we are using the entire region of Eu and 
not removing the resonance by the cut Ej. > E^q as in the case of cut 1 . Consequently a broader resonance 
and hence a larger value of \h\ gives a larger value of D and vice versa. 

Figure 4 shows the squared amplitude for B° -»• D*+D*-R s and B° -»• D* + D*~K S as a function of 
the variable y for h = —0.4. As mentioned above the nature of the two curves reflects the fact that in the 
process B° — ► D* + D* s '{~ — > D* + D*~K° the kaon emerges most of the time closer to £>*~ than the D* + 
while the situation is reversed for B° decays. 

Fig. 5 shows the plot of the function Gq,G c , G s \ and G S 2 as a function of y for = 0.6 GeV and for 
h = —0.4. From the figure we see that the functions Go and G s \ are symmetric in y while G c and G S 2 are 
antisymmetric in y. This follows from the absence of direct CP violation as shown in Eq. (19-22). 

In Fig. 6 we show the the decay distribution dT/dE^ versus the kaon energy E^. For small values of 
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Figure 5: The function Go,G c , G s \ and G S 2 as a function of y for = 0.6 GeV and h = —0.4 

Ek the decay distribution shows a clear resonant structure which comes from the pole contribution to a\ res 
with the excited J p = 1 + intermediate state. Therefore, examination of the D*K S mass spectrum may be 
the best experimental way to find the broad 1 + p-wave D s meson and as mentioned in the previous section 
a fit to the decay distribution will measure its mass and the coupling. 

In Fig. 7 we show the functions Go, G c , G s \ and G S 2 integrated over the y > as a function of h. Jq 
and J c refer to the integrated Go and G c functions while J sl and J s2 refer to the integrated G s \ and G S 2 
functions. As already mentioned, restricting the integration range to y > allows a probe of the cos (2/3) 
term in the time dependent rate for B°(t) — ► D^* +S) D^*^K S decays. It is clear from the figure that a 
broader resonance is more favorable to probe G S 2 which is the coefficient of the cos(2/3) term. 

In summary, we have studied the possibility of extracting sin(2/3) and cos(2/3) from time dependent 
B° —>■ D^D^ K s decays. These decays are expected to have less penguin contamination and much larger 
branching fractions than the two body modes B° — > D^D^ . Using HHCHPT we have calculated 
the branching fractions and the various coefficient functions that appear in the time dependent rate for 
B° — > £)(*+) D(*~^ K s . We also showed that a examination of the D*K S mass spectrum may be the best 
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Figure 6: The decay distribution dT/dEk versus the kaon energy E^. 
experimental way to find the broad 1 + p-wave D s meson and measure its mass and coupling. 
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Figure 7: The functions Gq,G c , G s \ and G S 2 integrated over the y > as a function of h. Jo and J c refer 
to the integrated Go and G c functions while J s ± and refer to the integrated G s \ and G S 2 functions. 
The values of the integral can be obtained by multiplying by Tb where is the width of the B. 

Appendix A 



In the Standard Model (SM) the amplitudes for B — > D^D^K S are generated by the following effective 



Hamiltonian |24], [2§] : 



H eff 



Gi- 



lo 



-±[v fb v f * q ( Cl ot f + C2 o q 2f ) - £(K b KVr + v cb v: q 4 + V tb V t * q ct)Of] + B.C. , (A. 1) 

where the superscript u, c, t indicates the internal quark, / can be u or c quark, q can be either a d or a s 
quark depending on whether the decay is a AS = or AS = —1 process. The operators 0\ are defined as 



Of, = q a ^Lfpf^Lb a , O g 2f = q ltl Lff^Lb, 



u 7,9 



q-y^Lbq'-i^L(R)q' , 0| 6 = q a ~i ^Lb p q'^ ^L(R)q' a , 
3 3 

-Ti^Lbe^i l^R{L)q' , O§ il0 = -g a 7 At L6 (3 e (? /^7 Ai i?(L)^ , 



(A. 2) 



where -R(-L) = 1 =t 75, and </ is summed over all flavors except t. Oif^f are the current-current operators 
that represent tree level processes. O3-Q are the strong gluon induced penguin operators, and operators 
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Oj-io are due to 7 and Z exchange (electroweak penguins), and "box" diagrams at loop level. The Wilson 
coefficients c{ are defined at the scale /i ~ m& and have been evaluated to next-to-leading order in QCD. 



The q are the regularization scheme independent values obtained in Ref. ||26[|. We give the non-zero c; 
below for m( = 176 GeV, a s (mz) = 0.117, and /i = nif, = 5 GeV, 

Cl = -0.307 , c 2 = 1.147 , 4 = 0.017 , c\ = -0.037 , 4 = °- 010 . 4 = -0.045 , 

4 = -1.24 x io~ 5 , 4 = 3 - 77 x 10 ~ 4 > 4 = -0.010 , 4o = 2 - 06 x 10-3 > 

4,'5 = -cS/JV e = W, <' 9 c = Pr, 4^ = (A. 3) 

where iV c is the number of color. The leading contributions to P\ e are given by: P\ = (§^) c 2(ip + 
G(mi,fM,q 2 )) and = (-^)(iV c ci + c 2 )(^ + G(mi, fj,,q 2 )). The function G(m,fj,,q 2 ) is given by 

Gfm, /x, g 2 ) = 4 /* gfl - jOln "** ~ ^ ~ ^ dx . (A. 4) 

Jo m 

All the above coefficients are obtained up to one loop order in electroweak interactions. The momentum 
q is the momentum carried by the virtual gluon in the penguin diagram. When q 2 > 4m 2 , G(m, fj,,q 2 ) 
becomes imaginary. In our calculation, we use m u = 5 MeV, = 7 MeV, m s = 200 MeV, m c = 1.35 



GeV f27|, |3|. 



In the factorization assumption the amplitude for B — > D^D^ K s can now be written as 

M = Mi + M 2 + M 3 + M 4 (A. 5) 



where 



Mi = -^Xi < J D ( * ) ^ s |s 7 ^(l-7 5 )c|0>< J D w |c7 Ai (l-7 5 )&| J B> 
a/2 

M 2 = ^fx 2 <D { * ) D^\c^{l-^)c\Q><K s \s 1 ^l-^)b\B> 
v2 

M 3 = ^X 3 <D { * ) D^\c^{l + -f 5 )c\0><K s \s-f^l--f 5 )b\B> 
V2 

M 4 = ^X 4 <7J ( * ) ^|s(l + 7 5 )c|0>< J D ( * ) |c(l-7 5 )6| J B> (A. 6) 

v2 

(A. 7) 
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where 



X 2 = V c (d + §-)+B 3 + ^-B A + B 9 + -L^io 



-^3 = B$ + — Bq + B-j + —B$ 




X A = -2 (^-B 5 + B 6 + ^-B 7 + B 8 



(A. 8) 



We have defined 



Bi = - E ^ 



(A. 9) 



q=u,c,t 



with 



(A. 10) 



In the above equations iV c represents the number of colors. It is usually the practice in the study of 
two body non-leptonic decays to include non-factorizable effects by the replacement N c — > N e ff. Since it 
is not obvious that N e ff for two body non-leptonic decays is the same for non-leptonic three body decays 
we will use iV c = 3 in our calculation. 

As already mentioned, we expect the contribution from penguin diagrams to be small and so as a first 
approximation we will neglect M3 and M4. Furthermore, from the values of the Wilson coefficients 
given above in the previous section it is clear that the amplitude M2 is suppressed with respect to Mi with 
the Wilson coefficients associated with M2 being about 7 % of the Wilson coefficients associated with M\. 
We also note that the currents < D^K S \ 57^ (1 — 7 s ) c | > and < K$\ s 7 M (1 — 7 s ) b \ B >, which appear 
in Mi and M2 respectively, receive contributions from both the contact terms and the pole terms. For the 
former current the pole terms are proportional to \/{Ek — Sm) while for the latter the pole term goes as 
\j(EiK + Sm). This also leads to a further suppression of M2 relative to M\. We therefore neglect M2 and 
only retain Mi in our calculation. We will also neglect CP violation in the K° — K° system and so (with 
an appropriate choice of phase convention) we can write 



K s = 



K°-K° 
V2 
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To calculate the various matrix elements in M\ above we use Heavy Hadron Chiral Perturbation Theory 
(HHCHPT). In HHCHPT, the ground state (j p = ±~) heavy mesons are described by the 4x4 Dirac 



matrix 



(A. 11) 



where v is the heavy meson velocity, P* M and P a are annihilation operators of the 1 and Qq a mesons 
(a = 1, 2, 3 for u, d and s): for charm, they are D* and D respectively. The field H a is defined by 



Ha = 7°# V 



Similarly, the positive parity 1 + and + states (j p = ^ + ) are described by 

(1 + 



S a 



(A. 12) 



In the above equations v generically represents the heavy meson four-velocity and D*^ and D are 
annihilation operators normalized as follows: 



<0|D|cg((T)) = 7^ 
<0|Z>^|cg(l-)) = e^^. 



(A. 13) 



Similar equations hold for the positive parity states and Dq . The vector states in the multiplet satisfy 
the transversality conditions 



v »d; = v^D*^ = 0. 



For the octet of the pseudo Goldstone bosons, one uses the exponential form: 



where 



and f w = 132 MeV. 



iM 

S = exp — 

JTC 



( Jl*° 



M = 



r + 



7T 

K~ 



7T 



o 



V J 



(A. 14) 



(A. 15) 
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The lagrangian describing the fields H , S and £ and their interactions, under the hypothesis of chiral 
and spin- flavor symmetry and at the lowest order in light mesons derivatives, is|14}|: 

f 2 

J 7T / 



-Tr 



+ iH b V»D^ ba H a 

+ Tr [S b (i vPD^ - 5 ba A)S a ] + igTr [H hW A^H a ] 
+ ig'Tr[S b i^ 5 A^S a ] + ihTr[S b ^ 5 A% a H a ] + h.c. 



(A. 16) 



where Tr means the trace, and 



D^ha = had,, + Vf.ha = 5 ba d^ + - (tfdpt; + £8^ 



ba 



(A. 17) 
(A. 18) 
(A. 19) 

£ = £ 2 and A is the mass splitting of the S a states from the ground state H a . 

The currents involving the heavy b and c quarks, Jy =< D* (ei,pi)\c^^(l—^)b\B(p) > can be expressed 



in general in terms of form factors [28] 
-2iV(q 2 



J V 



m + mi 



£^ a0 e* luPaPl p - (m + m 1 )A 1 {q 2 )e 1 x + ^^-e\ ■ q(p + Pl f 



m + mi 



+ 2m 1 A 3 (q 2 )X^q»-2m 1 A ( q 2 )X^q fl 



(A. 20) 



with 



A 3 (q 2 ) 



m + mi 
2ni\ 

4,(0) = 4,(0) 



A\{q 2 



2ni\ 



-A 2 (q 2 ) 



(A. 21) 



where q = p — p\ is the momentum transfer and m and ni\ are the masses of B and D* . In the heavy 
quark limit the various form factors are related to a universal Isgur-Wise function £(y • v{) where v and v\ 
are the four velocities of the B and the D* meson. One can write 



myjni\^(v ■ Vi) 



-is 



liva/3 * 



e*vV a vi/3 + e* • v - ef{v ■ vi + 1) 



(A. 22) 
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The weak current = g a 7 M (l — 7s)Q can be written in the effective theory as 

>(1 -75)^ fe + J (A. 23) 



a 2 



where Jq is the heavy meson decay constant. One can therefore write 

<D*(e 2 ,v 2 )K°\s^(l- l5 )c\0> = i m2f ^ (A. 24) 

JK 

Appendix B 

The total amplitude for B°(v,m) — > Z)* + (e 1 , v + , mi)D*~ (e 2 , w_ , mi)K s (pk) can be written as 

a = a non - res [l — P\] (B. 1) 

and the total amplitude for B°(v,m) — > Z?* + (ei,w + ,mi)D*~(e2,w-,mi)i^ s (pfc) can be written as 

a = Orion-res[l - -P2] ( B -2) 

with 



/d*{ / m*' /ipfc • w_ 



f D . V mi ^ . y _ + mi _ m « + l_D£j 



(B. 3) 



P 2 = M/Z ^± __ (B .4) 

/d* V m i (p k . v+ + mi _ m *> + 1^1) 



2 

In the boosted frame we can write 

_ K^gLv 

m\ 
m\ 

where E' k and p' k are the energy and the magnitude of the momentum of the kaon in the boosted frame, 
E'± and p'± are the energies and the magnitude of the momenta of the D*^ 1 in the boosted frame and mi 
is the D* mass. In the boosted frame we have the following relations 

E' k = 7 (E k -P-p- k ) (B. 7) 
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1 E l~ m l 

1 ET- 



2 

( ■, , / / / 

E k + 



El - vn 



m 2 (l -) 2 



(B. 



^ = V^ 2 - ml (B. 9) 



-/_ = - ml (B. 10) 



Pfc = Pb 
P+ = V 

E ' + = E'_ = E ' B ~ E ' k (B. 11) 

where and pt are the energy and magnitude of the momentum of the K s is the B rest frame, E' B and 
p' B are the energy and magnitude of the momentum of the B in the boosted frame and m, mi and mt are 
the B, D* and K s masses. 

Note from the above relations that Pi and P2 can be expressed in terms of E^ and y and P\(y,Ek) = 
Pb(-y,E k ). 

Squaring the amplitudes and summing over polarizations one can write 

|fl| 2 = |araora-res| 2 |l — -Pl| 2 (B. 12) 

|a| 2 = |a rao ?i-res| 2 |l — -P2I 2 (B. 13) 

a*a = a* non - res a non - res (l — P2)* (I — P\) (B. 14) 

where 



\ a non— res | 



: -' = k 2 [-x 2 + 2(2xix 2 + x 2 )x + 2x 2 - x\ + 4xi + 2] (B. 15) 



|2 _ ,„2r „2 , n ^o„, ™ i „ \„ i o™2 „2 
Inon— res | 



non— res^non— res 



rc^[-x 2 + 2(2xix 2 + xi)x + 2x^ -x( + 4x 2 + 2] (B. 16) 

= K 2 [x 2 + (xi + x 2 -2)x + 2xi + 5xix 2 + 2x 2 + 2 + 0(^/m 2 )] (B. 17) 



where 



k = —j=V c [ + c 2 ) \fm^fm{m\ 



V2 

X\ = v ■ v + 



(t +c2 ) ^ 

E' B E' + -p' B p' + y 

e' b e'_ + p' B p'-y 



x 2 = v ■ V- = 

mmi 



25 



E' + E'_ +p'p'_ 
X = v + ■ V - = — 



The double differential decay distribution for the time independent process 

B°(v,m) -> D* + (e 1 ,v + ,m 1 )D*-(e 2 ,v_,mi)K s (p k ) 

can be written as 

1 dT f(y,E k ) 
TdydE k Ifi y jEk) Ej^ dydEk 

f(y,E k ) = [-x 2 + 2(2xiX2 + x 2 )x + 2x\-xl + 4xi + 2]|1-Pi| 2 (B. 19) 

where p^. and are the magnitudes of the three momentum of the K s and D* + in the boosted frame. 
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